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Fermi arcs and isotope effect of the magnetic penetration depth 
in underdoped cuprates 
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Abstract. - The isotope coefficient ft of the magnetic penetration depth in the superconducting 
state is studied at T = for a d-CDW and a nodal metal model. Disregarding superconductivity 
the Fermi surface of the first model possesses arcs whereas the second model has no arcs. We show 
that a large increase of ft in the pseudogap region is generically incompatible with Fermi arcs in the 
pseudogap state. Thus only the second model shows a large increase of ft with decreasing doping. 
The required electron-phonon coupling is small and compatible with first-principles calculations 
based on the local density approximation (LDA). 



Underdoped high-T c cuprates show properties which are 
not expected to occur in Fermi liquids. pQ There is an en- 
ergy gaplike feature (pseudogap) seen already well above 
the transition temperature T c which increases with de- 
creasing doping. At the same time the large Fermi surface 
of the normal state at large dopings transforms in the un- 
derdoped region into arcs around the nodal direction and 
becomes gapped near the antinodal points. |2J The tem- 
perature dependence of the length of these Fermi arcs is 
presently controversely discussed. In particular, it is de- 
bated whether for T — > and in the absence of supercon- 
ductivity their lengths would approach zero [3] or a finite 
value. [1] 

Another unexpected feature of underdoped cuprates is 
the observed large isotope effect ft for the magnetic pen- 
etration depth, ft is very small in the overdoped region, 
increases strongly with decreasing doping, and may then 
reach values of the order of one. [51[S] Such large values 
cannot be explained by Eliashberg theory where ft is prac- 
tically zero in agreement with a recent experiment in the 
strong-coupling superconductor MgB 2 [7]. 

The aim of this Letter is to show the strong interrelation 
between the pseudogap and Fermi arcs at T=0 and ft in 
the underdoped regime. Anomalously large isotope coeffi- 
cients due to the presence of the pseudogap have been sug- 



gested previously. [8] For the coefficient a, describing the 
change in T c due to isotope substitution, detailed quan- 
titative calculations [HEl] yielded large enhancements of 
a even for a small electron-phonon (EP) coupling due to 
the presence of a pseudogap in rough agreement with the 
experiments. Similar calculations for ft, presented in this 
Letter, yield a more complex picture. An enhancement of 
ft does not only require a pseudogap but it is generally in- 
compatible with Fermi arcs at T — formed by infinitely 
sharp quasiparticles. To show this we first study a d-CDW 
model which can be derived from the t-J model within 
mean field theory. Such a model is generic for models 
where the pseudogap is associated with long-range order 
in the particle-hole channel yielding Fermi arcs at T = 
due to imperfect nesting. We show that this model leads 
to neglegible values for ft in the underdoped region be- 
cause some integrals perpendicular to the arcs diverge if 
the superconducting gap tends to zero. One way to avoid 
these divergencies is provided by the nodal metal model 
which has no arcs. Explicit calculations show in this case 
indeed a large increases of ft and a ratio of ft /a of about 2 
in the underdoped region in agreement with experiment. 

Our calculation is based on the large N limit of the t — J 
model (N is the number of spin components) [llj which 
represents a model with competing rf-CDW and d-w&ve su- 
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perconducting order parameters. [T^] For calculating /3 we 
add as in Ref. [TU] a phonon-induced interaction V(k— k') 
between electrons. In the following only the d-wave part of 

V is required which is obtained by replacing V(k— k') by 
4V7(k)7(k') which defines the d-wave coupling constant 

V for a phonon-induced nearest neighbor interaction. 7(k) 
is equal to (cos(fc x ) — cos(fc y ))/2. 

The isotope coefficient j3 for the phase stiffness A is de- 
fined by f3 = | gfa^ ■ A- i s related to the magnetic pene- 
tration depth A by A = c 2 /(47re 2 A 2 ) where c and e are the 
velocity of light and the electronic charge, respectively. 
lod is the Debye frequency which is assumed to be pro- 
portional to M -0 - 5 where M is the ionic mass. At zero 
temperature only the diamagnetic term contributes to A 
in the superconducting state which is given by 
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Here and in the following we use a reduced zone scheme an- 
ticipating that the original high-temperature primitive cell 
doubles for a d-CDW ground state. Thus the momentum 
k in Eq. ([T]) runs only over half of the high-temperature 
Brillouin zone whereas the index a = 1, 2 counts the orig- 
inal and the backfolded electronic branches. According to 
Peierl's substitution rule X<*(k) are the electronic eigen- 
states in the normal state renormalized by many-body in- 
teraction within mean-field theory. Thus they may de- 
scribe a pseudogap in the particle-hole channel but are 
unaffected by the superconducting gap. In our case they 
are given by 
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Xl , 2 (k) = ^±^ !ki 



4$ 2 (k), 



(2) 



with e±(k) = e(k)±e(k— Q), where Q = (ir, n) is the wave 
vector of the d-CDW. $(k) is the amplitude of the d-CDW 
and equal to $o7(k) with 7(k) = (cos(fc a ) — cos(fc y ))/2. 
The bare electronic energies e(k) are those of the t-J 
model in the large N limit counted from the chemical po- 
tential jit, i.e., e(k) = — 2(St + rJ)(cos(k x ) + cos(k y )) — 
4t'S cos(k x ) cos(fcy) - /i, with r = 1/N C ^ q cos(g x )/(e(q)). 
/ is the Fermi function, 6 the doping away from half- 
filling, J the Heisenberg coupling constant and t and t' 
are hopping amplitudes between nearest and next nearest 
neighbors on a square lattice, respectively. N c is the orig- 
inal number of primitive cells. n Q (k) denotes the electron 
density with momentum k and is given by 



n a (k) = 2Tj2G[ a 1 \k 7 iu Jn )e 



(3) 



where the prefactor 2 accounts for the spin degeneracy, 
rj is an infinitesimally small positive quantity and u n the 
Matusbara frequency (2n+l)irT. is the element (1,1) 



of the Green's function matrix G^ Q ^(k, iuj n 
inverse as, 



defined by its 
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Fig. 1: D-wave projected density of states iVd(w) in the pseudo- 
gap state without superconductivity for the the d-CDW (upper 
diagram) and the nodal metal model (lower diagram) at T — 0. 
Left and right insets show Fermi lines and Nd(0), respectively, 
at different dopings. 



The superconducting gap function A(k, n) consists of a 
contribution due to the <i-wave part of the Heisenberg and 
one due to the <i-wave part of the EP interaction, 

A(k,n) = -fVjAi + VFe(w D -K|)A3)>^y(k), (5) 



where is the step function. In Eqs. d}-© we used 
the fact that the phonon-induced reduction of the quasi- 
particle weight is independent of the ionic mass |10j and 
thus may be neglected in the following. 

The self-consistency equation for A splits up into two 
equations for Ai and A3, reading 



(1 + Fu)Ai + F 12 A 3 = 0, 
Fi 2 Ai + (l + F 22 )A 3 =0, 



(6) 
(7) 



with 
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(9) 



F22 is given by the expression for F±2 if one replaces there 
VJVby V. 

One important ingredient of the theory is the d-wave 
projected density Nd(uj) of the d-CDW model. It is shown 
in the upper diagram of Fig. [T] together with Fermi lines 
in the left inset. For the calculation we took the values 
J = 0.3 and t' = —0.35 measuring all energies in units of 
t. The right inset in Fig. Q] shows that A^(0) is nonzero 
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doping 



Fig. 2: (Color online) Isotope coefficients a, ao, ft, and /3j no for 
the d-CDW model as a function of doping using V = 0.04 and 
ujd = 0.16. Inset: Suppression of T c as a function of doping 
for the d-CDW model due to the evolving pseudogap. [TO] 



at finite 8 due to finite arcs and approaches zero only in 
the limit 8 — > 0. 

Recent LDA calculations [131 E] showed that the EP 
interation and in particular its d-wave part [T5] is very 
small in cuprates. This allows to simplify Eqs.© and © 
by keeping only the leading terms in V. In a first step one 
considers V = so that A3 = and determines in the d- 
CDW case self-consistently Ai = Ai and <I>(k) similar as 
in Ref. [TD]. Writing then Ai = Ai + A2 one may neglect 
F22 in Eq.0 and obtains from this equation 



A 3 = -F 12 (Ai ) 0,w rj )A 1 . 
Eq.(0) yields, again in leading order in V, 

aFii(Ai,0,w D ) 



(10) 



dAi 



(11) 



After taking the derivatives in Eq. (fTT]) Ai and A3 should 
be put to Ai and 0, respectively. Assuming that A de- 
pends on Ud only via Ai and A3 and noting that Eqs. 
(TTUj) and (HU hold for a general ujo one may easily form 
the derivatives of A2 and A3 with respect town and thus 
obtain j3. 

In principle there are two more sources for a depen- 
dence of A on M: The bare one- particle energies e(k) may 
depend on M. Such a case may occur if the EP inter- 
action can no longer be described by Eliashberg theory 
because strong polaronic or non-adiabatic effects play an 
important role. However, there is overwhelming evidence, 
both from experiment [16] and theory [13j[14] , that this is 
not the case. For instance, recent angle resolved photoe- 
mission (ARPES) experiments on nearly optimally doped 
Bi2212 [TDJ show that the band width does not change 
within the experimental error when substituting O 16 by 
O 18 . Using the polaronic model [T7] an observed value of 



(3 ~ 1 would require a 12 per cent change of the electronic 
band width which is not observed. Another possibility 
for a M dependence of A could be due to the d-CDW gap 
$(k). Recently we have shown [TUJ that the onset temper- 
ature T* for the <i-CDW gap shows practically no isotope 
effect which should also hold for the T — d-CDW gap 
$(k). 

Numerical results for the isotope coefficients a, ao, and 
/3, which are related to T c , the T = superconducting gap 
and the phase stiffness, respectively, are shown in Fig. [2] 
for the d-CDW model. The calculation of a has been 
described previously [10) . ao is defined by 



a 



ujd dAi 
2A~i du D ' 



(12) 



Outside of the pseudogap state, i.e., for 8 > S c ~ 0.145, all 
three coefficients are very small. Here a and ao are much 
smaller than the BCS value of 1/2 because superconduc- 
tivity is mainly determined by J and not by V because of 
the small employed value for V. Entering the pseudogap 
state a and ao nevertheless increase strongly with decreas- 
ing doping due to the evolving pseudogap as explained in 
detail in Ref. [10]. In contrast to that (3 does not change 
much below 8 C and remains practically zero throughout 
the underdoped region. Considering only the underdoped 
region, i.e., 8 < 5 C , one could also use the ratio T c /T c0 as 
an independent variable instead of 8. The curves in Fig. 
[2] would not look much different in such a plot because T c 
depends rather linearly on 8 in that region, as shown in 
the inset of Fig. [2] 

To understand the surprising behavior of /3 we present 
now approximate analytic expressions for /? and various 
other quantities which become asymptotically exact at 
small values of the superconducting gap but are also rather 
accurate over the whole doping region. Evaluating Eq.(UJ 
for an infinite cutoff and neglecting momentum derivatives 
of the order parameter we find 
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The main contribution in the sums over momentum in 
Eqs. ([TBI and (fT4"]l comes from the region near the arcs 
which consists of all points {ki? a } satisfying %(kjr a ) = 0. 
Correspondingly, we split these sums into a part parallel 
and a part perpendicular to the arcs. Assuming that the 
electron dispersion perpendicular to the arcs can be ap- 
proximated linearly we find that the integrations perpen- 
dicular to the arcs diverge for small A(k) cancelling, for 
instance, the prefactor A 2 (k) in Eq. (fT3]) . The remaining 
integration parallel to the arcs in A and X can be written 
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as surface integrals along the arcs yielding 

A^/jS^VU^))'. (15) 

Ot 

and X — A/ 3, or equivalently, 

P/ao = 0. (16) 

Eq. fToT) explains the tiny values for f3 in Fig. [5] obtained 
by numerical evaluation of the exact expressions for /3. 
The small values for /3 throughout the underdoped region 
reflect directly the divergencies encountered in the mo- 
mentum integrations perpendicular to the arcs. If there 
would be no divergencies the right-hand side of Eq. (fT4|) 
would smoothly approach the right-hand side of Eq. (fl"3|) 
for A(k) — > 0. Thus we would obtain (3/ao — 2 and an 
increasing /3 towards small dopings both of which would 
be in excellent agreement with experiment. This suggests 
that the above mentionned divergencies actually do not 
occur in the high-T c superconductors so far studied with 
respect to j3. 

One model which has a pseudogap but no arcs and thus 
no divergencies in momentum sums is the nodal metal 
(NM) model [3j[9]. Its dispersion is given by 

Xl , 2 (k) = ±Ve 2 (k) + $ 2 (k), (17) 

which formally can be obtained from Eq. ^ by putting 
e+(k) = 0. The NM model does not allow to determine 
$o but considers it as a parameter. Following Ref. [9] it is 
convenient not to use $o as an independent variable but 
the ratio T c /T c0 which describes the reduction of T c due 
to ^(k) relative to the maximum transition temperature 
T c o at 6 = S c . Dependencies on T c /T c q may be interpreted 
as doping dependencies due to the monotonic relation be- 
tween T c and S in the underdoped region, sec the inset of 
Fig. [2] in case of the d-CDW system. The lower diagram 
in Fig. Q] shows N d (u) for the NM model for 5 = 0.10 and 
two values for T c /T c q. Nd(0) is zero in all cases because 
of the absence of arcs. 

The calculation of (3 proceeds in the same way as for the 
d-CDW model. Fig. [3]shows numerical results for a, a , /3, 
and the ratio /3/ao using V = 0.04 and lo d = 0.16. Un- 
like in Fig. [5] (3 increases strongly with decreasing T c /T c q 
and the ratio /3/ao is near two over a large part of the 
underdoped region, approaching exactly 2 for T c /T c q — > 0. 
Since a and ao are of similar magnitude this implies also 
PI a ~ 2 which is exactly the value found in experiment. 
The diagram also shows that /3 assumes values of 1 and 
larger already for moderate reductions in T c . This is rather 
astonishing because it means that the small EP coupling 
constant found in LDA calculations [T3HT5] is able to pro- 
duce the large value for j3 of Fig. [3j Varying the Debye 
frequency over a large frequency region does not change 
much the curves in Fig. [5] and Fig. [3] Similar as in the 
case of a discussed in Ref. [10] ao and /3 do not show 
anomalies if the phonon frequency and the pseudogap are 
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Fig. 3: (Color online) Isotope coefficients a, ao, /3 and the ratio 
/3/ao for the nodal metal model as a function of T c /T c o for a 
fixed doping 8 — 0.10 using V = 0.04 and uid = 0.16. 

close to each other. Actually, for all phonon frequencies 
smaller than about 0.16 the pseudogap crosses the phonon 
at some doping without any effect on the curves showing 
that the large increases in the underdoped region are not 
due to resonance effects between the phonon and the pseu- 
dogap. 

Important details of Fig. [3] such as the behavior of the 
ratio /3/ao can be understood again by an approximate 
analytic evaluation of the sums over momenta. Since the 
NM model has no arcs one may expand e(k) around the 
high-temperature Fermi line given by e(k) = 0. Splitting 
again the momenta sums into perpendicular and parallel 
parts, taking $(k) right on the Fermi line and linearising 
e(k) the integration perpendicular to the Fermi line can be 
carried out and is always finite. As a result (3 and X can 
again be written as surface integrals. Due to the absence 
of divergencies X — >• A for A (k) — » yielding /3/ao — > 2. 
This result is rather general and just expresses the absence 
of divergencies. 

Many features of the NM model agree with the experi- 
ment, for instance, the observed ratio /3/a ~ 2 over most 
of the underdoped region and the large increase of a and 
j3 with decreasing T c . Obviously it has also short-comings. 
For instance, the lower diagram of Fig. 1 implies a pseu- 
dogap of only about 0.02 or 10 meV for a decrease T c /T co 
of 0.3 which is unrealistically small. The more fundamen- 
tal question, however, is whether this model with perfect 
nesting in the particle-hole channel is applicable to the 
cuprates at all. ARPES experiments in extremly under- 
doped Bi2212 [TB] found an excitation spectrum of the 
pseudogap phase without superconductivity which is very 
similar to the NM model. These data, however, have been 
interpreted in terms of a superconducting order parameter 
with finite phase correlation length. Since superconduc- 
tivity does not renormalize the one-electron spectrum en- 
tering the phase stiffness the resulting /3 would be equally 
small as in a superconductor without pseudogap and thus 
disagree with experiment. Nodal metal behavior can, how- 
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ever, be realized to some extent for an incommensurate d- 
CDW state. The nesting condition e+(k) = of the NM 
model can in this case fulfilled in limited regions around k 
points on Fermi lines which are connected by nesting vec- 
tors and exhibit parallel tangents to the Fermi line. As- 
suming four nesting vectors of the form Q = (tt, tt ± Si nc ) 
and (tt ± Si nc , tt), determining Si nc from the nesting con- 
dition and correlating only pairs of states with maximal 
nesting we obtained the dotted line f3i nc in Fig. 2 for /3. 
It shows the desired increase towards low doping but the 
calculated ratio j3i nc /ao,i nc (not shown in Fig. 2) deviates 
substantially from 2 due to divergencies at Fermi lines 
near the antinodal point. It also has been shown [19] that 
superconductivity (not taken into account above in de- 
termining Si nc ) suppresses Si nc strongly, suggesting that 
incommensurability is not causal for the large observed 
values for (3. 

The above results show that the presence of arcs (or, 
more general, a finite length of the Fermi line) at T = 
in the pseudogap state (i.e., without superconducitvity) 
are generically incompatible with large values of /? in the 
underdoped region. This holds for all models where the 
pseudogap phase has long-range order and arises in the 
particle-hole channel which necessarily leads to arcs in 
two-dimensional models due to imperfect nesting. The 
(i-CDW model is one example. It also holds for short- 
range models 20 for the pseudogap which are Fermi liq- 
uids, i.e., where the self-energy at T = can expanded in 
powers of the momentum and frequency. In this case the 
arcs are nothing else than the Fermi lines in the normal 
state without pseudogap and the arguments used above 
for Fermi arcs apply. Large values for /3 can therefore 
only expected if no arcs exist at all as in the NM model 
or that the divergencies in the integration perpendicular 
to the arcs in calculatling j3 are not present for some rea- 
son. One possibility could be a non-Fermi liquid ground 
state without infinitely sharp quasiparticles and a broad 
spectral function. 

In conclusion we have investigated the isotope coeffi- 
cient f3 of the magnetic penetration depth at T = using 
two models for the pseudogap. The calculation shows that 
ft depends severely on the presence or absence of Fermi 
arcs at T = in the absence of superconductivity. The 
experimentally observed large increase of (3 in the under- 
doped region is reproduced in the case of a nodal metal us- 
ing a small EP coupling consistent with first-principles cal- 
culations based on the local density approximation. The 
(i-CDW model is generic for models with two competing 
order parameters often used to interpret experimental re- 
sults. [2"T1l2"2"] The resulting /?, however, is small throughout 
the underdoped region and disagrees with experiment even 
if the EP coupling constant is increased by orders of mag- 
nitude. Responsible for the absence of an enhancement of 
/3 is the presence of Fermi arcs which lead to divergent k 
integrals perpendicular to the arcs in calculating (3. We 
hope that the discovered link between arcs and the iso- 
tope coefficient f3 motivates more experimental work on 



the doping dependence of /3. 
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